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ABSTRACT 

We explore the possibility that Milgrom’s Modified Newtonian Dynamics (MOND) is a mani¬ 
festation of the modification of inertia at small accelerations. Consistent with the Tully-Fisher 
relation, dynamics in the small acceleration domain may originate from a quartic (cubic) 
velocity-dependence of energy (momentum) whereas gravitational potentials remain linear 
with respect to mass. The natural framework for this interpretation is Finsler geometry. The 
simplest static isotropic Finsler metric of a gravitating mass that incorporates the Tully-Fisher 
relation at small acceleration is associated with a spacetime interval that is either a homoge¬ 
neous quartic root of polynomials of local displacements or a simple root of a rational fraction 
thereof. We determine the weak field gravitational equation and find that Finsler spacetimes 
that produce a Tully-Fisher relation require that the gravitational potential be modified. For an 
isolated mass, Newton’s potential Mr -1 is replaced by Mao log(r/ro) where ao is MOND’s 
acceleration scale and r o is a yet undetermined distance scale. Orbital energy is linear with 
respect to mass but angular momentum is proportional to AT 5 / 4 . Asymptotic light deflec¬ 
tion resulting from time curvature is similar to that of a singular isothermal sphere implying 
that space curvature must be the main source of deflection in static Finsler spacetimes possi¬ 
bly through the presence of the distance scale ro that appears in the asymptotic form of the 
gravitational potential. The quartic nature of the Finsler metric hints at the existence of an 
underlying area-metric that describes the effective structure of spacetime. 

Key words: galaxies: general - galaxies: kinematics and dynamics - dark matter 


1 INTRODUCTION 


Milgrom’s theory of Mod ified N ewtonian Dynamics (MOND) was put forward to model phenomenologically the velocity excess of galactic 
rotation curves (Milgrom 19830 . In contrast to the dark matter theory. MOND assumes that only luminous galactic matter is responsible 
for the observed velocities but that Newtonian dynamics breaks down where stellar motion involves accelerations smaller than a constant, 
ao, that is of order cH where c is the light speed and H the Hubble constant. Newton’s second law then assumes the modified quadratic 
form a 2 a^ 1 oc F™ where a and are the inertial and Newton’s gravitational accelerations respectively. As a phenomenological theory, 
MOND’s ability to fit galaxy rotation curves is well established i McGaugi] 20081 : Milgrom 2012 ). This success is also paralleled b y the 
theory’s inability to model the dynamics of galaxy clusters without requiring an unseen m atter additive i Famaev & McGaughl 2012 ). and 
explain the morphology of galaxy cluster encounters ' Clowe et al]2006l : lBradac et alj2008 ) that seem to strengthen the case for dark matter. 
These failings are not necessarily intrinsic to MOND’s basic phenomenology as much as they are to its available formulations that may not 
describe the physical theory from which the scaling relations, motivated mainly by flat galactic rotation curves, originate. 

In search for the physical foundation of MOND, the modification of Newton’s second law has been viewed as an indication of a 
modification of gravity as well as an indication of a modification of inertia. Modifying gravity entails searching for theories where motion is 
described by a = F g and F g is no longer the Newtonian gravitational acceleration but is related to it by the scaling relatio n F a = ( apF D 1 / 2 
for small accelerations. The first such example is the widely used low-energy lagrangian-based theory of Bekenstein & Milgrom ( 19841) that 
modifies the gravitational Poisson equation nonlinearly. Various covariant theories for the modification of gravity along these lines have been 
formulated subsequently and allowed the application of the MOND paradigm to gravitational lensing and large scale structure formation. 
Modifying inertia amounts to explaining how for small accelerations, the inertial resistance of a body to an external force becomes a 2 a^ 4 
instead of the usual Newtonian linear scale-free behavior. Gravity in the weak field limit is treated just like any other force and given by the 


* E-mail: namouni@obs-nice.fr 


© 0000 RAS 
















2 Fathi Namouni 


standard Poisson equation. Attempts at deriving a formulation for the modification of inertia have not been successful because deriving a 
lagrangian based theory satisfying Galilean invariance requires the theory to be non-local' Milgromlll994 ). Together with the nonlinearity of 
MOND, this property makes finding a working formulation as a modification of inertia more difficult. More recently, it has been suggested 
that MOND phenomenology may be interpreted as the manifestation of spacetime scale invariance for small accelerations I Milgronj 20091 . 


2014). 


In this work, we seek a formulation of MOND as a modification of inertia in a different way. We base our search exclusively on 
the observations that MOND was initially conceived for and explains best, namely, flat galaxy rotation curves. We assume that Milgrom's 
original phenomenological scaling a 2 a ( ^ 1 oc Fg and the implied quadratic inertial resistance do not necessarily reflect a new universal 
dependence on the acceleration a in general but only the modification of inertia for nearly-circular orbits from v 2 r~ l at large accelerations to 
v 4 r~ 2 ag 1 at small accelerations where r and v are the circular radius and velocity respectively. Although this inertial scaling when equated 
with Newton’s gravitational acceleration GMr~ 2 gives the flat rotation curves as originally sought by Milgrom, we show in this work that 
a metric description of this inertial change requires the gravitational acceleration to be changed to GMa^r -1 where G is the gravitational 
constant. Although our starting point is the original MOND scaling, the approach we lay down in this work is based on metric spacetimes 
and at no point do we refer to the notion of force. The phenomenological scaling serves only as a guide for the kind of dynamical behavior 
we seek. 

An inertial scaling proportional to the fourth power of velocity is indicative of a possible Finsler structure of spacetime. In effect, in 
Finsler geometry, the proper time interval is no longer the square root of a quadratic polynomial of local displacements as in Riemann 
geometry. Instead, the spacetime interval may be a general homogeneous function of order 1 with respect to local displacements that can be 
chosen to produce a quartic dependence. In this context, our scope of implementing the Tully-Fisher relation for nearly circular orbits around 
a spherically symmetric mass distribution translates into searching for the simplest static isotropic Finsler metric that produces an inertia 
proportional to the fourth power of velocity for circular gravitational motion. 

What we seek to do for MOND in a Finsler spacetime is analogous to the following programme: knowing Kepler’s third law, one 
asks whether it is possible to write down the Riemann metric of a static isotropic mass distribution where inertia is proportional to the 
second power of velocity. If the mass distribution is concentrated in a single point, the sought-for result is the Schwarzschild metric. It is 
elementary, however, that Kepler's law does not define completely the Schwarzschild metric. Writing the Riemann spacetime interval as 
dr 2 = b(r)dt 2 — a(r)dr 2 — r 2 dcj> 2 in the standard spherical coordinates, one can derive the time potential b(r) from Newtonian theory but 
not the space potential, a(r) ( Gruber et al.ll988 ). The space potential may be determined from the behaviour of light deflection as both space 
and time curvatures contribute to the effect (equally). Incidentally, we remark that the time potential b(r ) obtained from Newtonian theory is 
exactly that which enters the Schwarzschild metric regardless of the assumed weak field approximation of Newtonian dynamics. 

Transposing this programme in the context of the MOND quartic scaling in a Finsler spacetime is more complex for two reasons. The 
first is the presence of more than two potentials in the metric’s expression (or more precisely in the Finsler function as explained in the 
next section). The second reason is the fact that light deflection is poorly known in the context of MOND as was mentioned above regarding 
gravitational lensing in merging galaxy clusters. We will however be able to determine the Finsler time potential that allows us to characterize 
MOND dynamics in the weak field approximation for massive motion. By considering light deflection in the same weak field approximation, 
we will be able to determine that asymptotic deflection by a static isotropic mass distribution is caused mainly by the space potential unlike 
that of general relativity. 

The paper is organized as follows: in section 2, we briefly review the physical applications of Finsler geometry and recall its main features 
that will be used in this work. In section 3, we derive the properties of Finsler spacetime metrics that produce MOND-like inertia, write down 
the corresponding weak field gravitational equation based on the Tully-Fisher relation at small acceleration and study its properties. Light 
deflection is analyzed in the weak field approximation in Section 4. We discuss our findings and their possible implications in Section 5. In 
the following, the gravitational constant and the light speed are set to unity. 


2 MODIFICATION OF INERTIA AND FINSLER SPACETIMES 


In classical physics, the non-relativistic linear dependence of momentum on velocity and the quadratic dependence of kinetic energy on 
velocity stem from the Riemannian structure of spacetime for which the proper time interval is given by the root of a quadratic form of 
local displacements as dr = ( g dx^dx 1 ') 1 ' 2 where g ^ is the metric tensor. That Nature has chosen this structure is obviously well 
established within the experimental range that current physics can reach. Nonetheless, different spacetime structures have been envisaged in 
various contexts: [Riemann 1 1868h considered the possibility that the spatial distance interval could be better described as a quartic polynomial 
of the local displacements rather than the standard quadratic form. The hypothesis was investigated extensively with proper time intervals 
of the for m dr = ( h Kl \ uv dx K dx x dx 1 1 dx l/ ) 1 ^ 4 , where h K \ ul/ is the gravitational 4-tensor, in the context of the metric tests of General 
Relativity (Rox burgh & Tavakol 19791 : Roxburgh 1991 . 1992al lhl). More recently, a non-Riemannian structure of spacetime appeared in an 
attempt to circumvent local Lorentz and CPT invariance by assuming that the local laws of physics are not invariant with respect to the 
full Lorentz group of transformations (Cohen & Glashow 20061) . Such limited invariance describes anisotropic spacetimes with the proper 
time interval dr = (ri u ,vdx a dx 1 ')( 1 - b ^ 2 (n 0 dx p ) b where |6| -C 1, is the Lorentz metric and the vector n p is position independent 
(Bogo slovsky & Goenneil 179991 : Gibbons et al.ll2007 ). String fluid motion in area-metric spacetimes that models localized massive particles 
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may also be described by a non-Riemannian quartic metric i Punzi et al]|2009al) . Lastly, exotic Finsler spacetimes have been considered as a 
possible premise for gravitation theories but their physical foundations are still lacking i Chang & Lil2008l:Er& Chan g||20 id). 

The relevant framework for describing such non-Riemannian modifications of inertia is Finsler geometry iBao et all 200ol) wherein the 
spacetime interval is written as dr = F(x, dx) and F is a positive homogeneous function of degree 1 with respect to local displacements dx 
in order to ensure proper time reparametrization invariance. A spacetime metric may be derived from the Finsler length function F as: 


JW( x,v) = 


d 2 F 2 


and f^(x,v)v IJ 'v 1 ' = F 2 (x,v), 


( 1 ) 


2dv^dv v 

where the latter property stems from the O-homogeneity of / M „ with respect to v. As the expressions of f are bound to be cumbersome 
because of the velocity dependence, it is common to use only the Finsler length function F. Geodesics of the Finsler spacetime are obtained 
from the usual length action: 


S = 




, dx\ dx M dx 1 ' 
\ X ’dXj~dX~dX 


dX. 


( 2 ) 


The existence of a well-defined causal structure in the Finsler spacetime requires the Finsler function to be of the form F = L 1 ^ n where L 
is a positive homogeneous function of degree n ^ 2 i Pfeifer & Wohlfarthll2011 : Pfeifeili2014 ). 

We approach the problem of the interpretation of MOND as a modification of inertia by considering a static spherical system in the small 
acceleration domain, and seeking the simplest static isotropic Finsler metrics that produce the Tully-Fisher relation, i> 4 oc M (r), for nearly 
circular orbits. It is clear from the outset that this condition will not define the Finsler metric entirely. The equations of motion we derive 
in the weak field approximation will be sufficient to ascertain the similarities and differences between the gravity and inertia modification 
interpretations of MOND in the small acceleration domain. In Section 4, we proceed further and examine the properties of asymptotic light 
deflection. 

On account of the causality requirement that the Finsler function must be of the form F = L 1 ^ 71 where L is n-homogeneous with n ^ 2, 
L must be a polynomial or a rational fraction of degree n with respect to local displacements. Limiting the search to isotropic spacetimes 
restricts n to even numbers. The simplest two possibilities are n = 2 and n = 4. For n = 2, L must be a rational fraction, the simplest 
being the ratio of two polynomials respectively of degree 4 and 2. Riemann geometry is obtained if L is a second degree polynomial. For 
n = 4, the simplest L is a polynomial of degree 4. On acco unt of spherical symmetry, F is function of only t he spherical radius r, dt, 
dr, and the solid angle element dfl = ( d6 2 + sin 2 ddtj/ 2 ) 1 ^ 2 i McCarthy & Rutz 1993 : Pfeifer & Wohlfarth 2012 ). Furthermore, in a static 
spherically symmetric spacetime, energy and angular momentum are constants of motion. This stems from of the Lagrange equations of the 
proper length action (0: d(d Vlx F)/dt — d x rF = 0. For isntance, as dtF = 0 then d Vt F = E is constant. 


3 STATIC ISOTROPIC FINSLER SPACETIMES FOR MOND 
3.1 Quartic metrics 


Quartic Finsler functions that produce MOND-like inertia must not have time and space cross terms. To see this, consider the gravitation-free 
F 4 = v 4 — atxV 2 v 2 — ait ) 4 where atx and a x are constants. The non-relativistic kinetic energy and linear momentum are expanded to 
fourth order in velocity as: Ek = d vt F — 1 = at x v 2 /A + (9aL/32 + 3a x /A)v 4 andpi = —d Vi F = atxVi/2 + (3a 2 x /8 + a x )v 2 Vi where 
i is a space index and Vi = dx 1 /dt. A quartic energy dependence requires atx = 0. Similarly, the most general static isotropic metric that 
modifies the small-velocity energy law from v 2 to v 4 is given by: 

4 


F 4 = 

r i 


dr V 

dA , 


= Pvt - (7 - a)vt - ( av 2 r + r 2 vfif, 


(3) 


where ty, = dx 11 /dA and the potentials a, /3 and 7 are functions of the radial variable r only. The coefficient of r 4 VQ was set to unity by 
choosing spherical coordinates that describe space as nested spheres of radius r for constant time t and radius r. We note that the original 
MOND scaling applies only to angular motion so that requiring the time and radius cross term to vanish is motivated by our view that inertia 
isotropy should take the same form in all space directions. 

To proceed further, we first consider the spacetime defined by F q with 7 = a 2 as the associated effective gravitational potential may be 
written explicitly. Hence: 

F q = \p 2 v t — av r - r vnj [P 2 v t + av r +r v Q J . (4) 

In order to apply the Tully-Fisher relation to circular orbits, we determine the timelike and null geodesic equations directly from the energy 
and angular momentum constants. Owing to spherical symmetry, orbital motion lies in a constant plane that we define as 6 = 7t/2. As a 
consequence, the solid angle fl is replaced by the longitude eft. Choosing the affine parameter A = r for massive particles, the constants of 
motion are given by E = d vt F q = pv/ and J = —d Vlj> F q = r 2 v^,(a u 2 + r 2 v |). Substituting the constants into the definition of proper 
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time 0 yields the orbital equations in terms of effective potentials as follows: 

14 1 J ^ 

avt = (B~' S E^ — l ) 5 -i-;- = — V e fr m for massive particles, 

r 2 (/3-3S3 - 1) 

2 _ 1 2 

av r = 8 e E 3 -—— = — Vefi o for photons. 

r 2 /3-3£3 


(5) 

( 6 ) 


Circular orbits for massive particles are the energy minima of the effective potential and satisfy: I4ff, m (fcirc) = 0 = 9rl4ft,m(fcirc). 
Combining the two equations and using the expression of the circular velocity v C irc(r) = rv we get: 


dr log/3(r) 


d Ccirr (r ) 

r [1 + Vcirc(r) 4 ]' 


(7) 


For comparison, the time potential, b(r), of a static spherically symmetric system with a Riemann metric dr 2 = bdt 2 — adr 2 — r 2 d(p 2 
satisfies the equation d r log b = 2 v 2 iic /r (1 + v 2 iic ). We may determine a simpler equation for the potential /3 by noting that the Tully- 
Fisher relation is verified only at low energy (v c i rc 1). Assuming further that spacetime is flat in the absence of matter, we may use in the 
weak field limit the fact that |/3 — 11 -C 1. Equation 0 then becomes d r p = 4 M (r)aor -1 where the Tully-Fisher relation, t; 4 irc = M(r)ao 
was used to eliminate the circular velocity. In terms of the matter density p = (Anr 2 )~ 1 d r M, the weak field gravitational equation reads: 


r 2 d r 



8ivp(r). 


(8) 


Following the same approximations but using the Kepler velocity relation instead of the Tully-Fisher relation, the time potential b[r) of the 
Riemann metric satisfies the Poisson equation r~ 2 d r (r 2 d r b) = 8np(r) in the Newtonian limit. To determine an expression of the time 
potential /3(r), regardless of the weak field approximation, requires the knowledge of the exact velocity relation or the full gravitational 
equation that replaces the Einstein equation in the small acceleration limit. Lacking both elements, we will continue our study of the Finsler 
metric using equation 0. The weak field gravitational equation should be able to yield the asymptotic behavior in the small acceleration limit 
(i.e. the dynamics far away from an isolated mass) in the same way the solution of the Poisson equation for the Riemann metric encapsulates 
the fundamental features of the dynamics of an isolated mass and incidentally yields the exact expression the time potential b(r) regardless 
of the weak field approximation. 

For a point mass, equation 0 yields /3(r) = 1 + AMao log(r/ro) where ?'o is a constant that may depend on M and ao. Whereas the 
constant r o does not influence qualitatively the properties of massive motion, its possible dependence on M and ao influences light deflection 
as we explain in Section 3. The energy and angular momentum of circular orbits are given as: 


- [1 - Ma a ~ 1 + M “» |3/4 + >’ 

J„, r . = (i _ Mo „ i Solog^A,)]»'- ~ P + ““I 1 ' 4 + ■ for Mao « 1. 

A first interesting result here is that modifying inertia through a non-Riemannian structure of spacetime requires a modification of 
gravity even though gravity remains linear with respect to mass at low energy. For comparison, the weak field modified gravity MOND 
potential is of the form ( Mao ) 1 ^ 2 log(r/ro), orbital energy is oc (Mao) 1 / 2 and angular momentum is oc (Mao^Vcirc. 

The timelike effective potential, Vefl.m, of a point mass is shown in Figure (l,left). Motion is always bound regardless of the unknown 

_ l 4 . 

potential a and the sign of orbital energy E — 1, as the condition /?“ 3 Es — 1 > Oimpliesr < ro exp[(i5—l)/Mao]. The potential’s details 
at small radii and in particular its unstable circular orbit that resembles that of the Schwarzschild timelike effective potential are unimportant at 
this stage as the small acceleration approximation is no longer valid and the metric structure should tend to that of a Schwarzschild spacetime. 
In this case, the condition that the time and space cross terms vanish in the Finsler function no longer holds if the Finsler spacetime is to 
describe all acceleration scales. 

The null effective potential, 14fr,o, of a point mass is shown in Figure (l,right). There are two equilibria, one stable and another 
unstable. Motion may be unbound. As with timelike motion, the unstable equilibrium is reminiscent of that of the Schwarzschild null 
potential. To ascertain the nature of the stable point we determine the radial motion extrema given by v r = 0. These occur at the distance 
r m = b 4 [l + AM log(r m /ro)] where M = Mao and b = J/E that may be termed the impact parameter. There are only two solutions that 
may be written as: 


Tm 


T m 


bM J 



4 -J_ \ 1 4 

rge m ^ 
b A M J 


bM i 




(9) 

(TO) 


where Wo and W -1 are the two branches of the real Lambert’s function whose argument satisfies —e 1 ^ — r^b 4 e ' 1 / M ^/M ^ 0 
implying b > ro The two solutions correspond to unbound motion, one of which 0 is unphysical since the photon is 
reflected at the larger r m value d. As the impact parameter increases, the unphysical solution tends to the zero of the time potential /3. 
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r Log 10 r 

Figure 1. Effective timelike (left) and null (right) potentials. For the fonner, a circular orbit (dashed) and an eccentric orbit (solid) are shown with 
(Mao, E, J) = (0.1.1.25, 0.05) and (0.1,1.132, 0.05) respectively. For the later, an unstable circular orbit (dashed) and an unbound orbit (solid) are 
shown with (Mag, E, J) = (0.1,17.15, 2.5) and (0.1, 9, 2.5) respectively. 


The apparent stable minimum is therefore never reached regardless of energy and angular momentum, implying that motion is unbound. We 
discuss light deflection in the next section. 

We now turn to the spacetime defined by the quartic length function F q with 7 ^ a 2 . The radial velocity is the solution of the polynomial 
equation of degree 6 in w 2 : 

1 


1 - 


(7 - a V 

p-iiEi -e 


2 

7W = 


( q 2 _ 7 )j 4 t[/ 3 3£3 - e - ( 7 -a 2 )v 4 ] 


r 4 (r^5 -e) 2 


/3-3E3 -e 


qJ 2 


r 2 (P 


( 11 ) 


where e = 0 or 1 for null or timelike geodesics respectively, and v r = dr/dr for massive particles. The new effective potential is an implicit 
function of the gravitational potentials. However, the weak field time potential /? is still given by equation ©. To see this explicitly, one 
may expand the energy constant of massive particles for both Vi -C 1 where Vi = dxi/dt, and for small mass by writing a = 1 + Ma, 
P = 1 + M/3 and 7 = 1 + M7, to find: 


E = 


3(v r + r vf) MP M p ^ ^ 4 a / A - 0 a\ 222 4 A-\ 

1 H-;-b -rrr + — [3(47 — 3 P)v r + 6(4a — 3 P)r v^Vr - 9 pr v ^ 


16 


16 


( 12 ) 


and J = r 2 t)^(t) 2 + r 2 -u|). weak field motion depends only on the time potential P through the third term in equation ( 1 1 2b except for 
the additional perturbation effects such as bound orbit precession given by the last term. Using the Tully-Fisher relation, the weak field 
gravitational equation ((8)1 may be recovered easily. The features examined earlier for the quartic metric with 7 = a 2 remain generally valid 
for massive particles when 7^a 2 . 


3.2 Rational metrics 

The general static isotropic rational length function may be written as: 


F 2 = 

V - 


p 0 v? - PiVr - p2r 4 Vft - p3r 2 vnvl + p A r 2 Vn'Vt. + Psvlvt 
v 2 + aiVr + a 2 r 2 t>o 


( 13 ) 


where v M = dx 11 /d\ and the potentials ap and pi are functions of the radial variable r only. Although unusual, a rational metric offers an 
interesting possibility with respect to quartic metrics. It may be possible to interpret the proper length, dr, as the ratio of quantity associated 
with a surface area to a length element. More specifically, the length function's numerator may be related indirectly to surface area through 
the Fresnel 4 -tensor of an area-metric spacetime (see Section 5 ). The obvious drawback of the rational metric is the larger number of unknown 
potentials that we reduce by requiring a quartic dependence for energy and cubic for angular momentum. 

As with the quartic metric, we set the orbital plane to 6 = 7 t/ 2 . The linear dependence of angular momentum J = —d v<j> F r on v$ is 
eliminated if ^4 = P00.2 as can be shown by simple algebra. Requiring that the quadratic radial velocity term vanishes (per velocity isotropy) 
in the energy constant E = d Vt F r leads to /3s = Poai. The length function then simplifies to: 

Pivt + «2r 4 t4 + p3r 2 vlv 2 


F 2 = p 0 Vt - 


v? + aitt 2 + a 2 r 2 tr? 


( 14 ) 


where P 2 = ot 2 was set by choosing the variable r to mean the physical spherical radius at constant radius and time as was done for the 
quartic metric. We note that the relations among ai, « 2 , Po, Pa and /3s no longer hold if the small acceleration domain is left and the Finsler 
function is to describe all acceleration scales and tend to a Schwarzschild metric. 
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The Tully-Fisher relation may be implemented by considering the weak field limit, v C irc <S 1, for which the constants of motion may 
be expanded to fourth order in velocity to yield the geodesic equations as follows: 

,.4 I a „2,.2„2 , „4„.4\ 

(15) 


p _ fli/2 , 3((3iVr + for 2 v?vl +a 2 r 4 vl) 

^ — Po + 


J = 


2/3 0 1/2 

^v^fovl + 2a 2 r 2 vl) 


Po 


1/2 


(16) 


where n = dx l /dt. Note that the potential ai is absent from the equations of motion. The solution for the radial velocity v r is unsurprisingly 
similar to that of the general quartic metric CD): 


2/3i 


1 - 


3 (fo - pl/Aa 2 ) v 2 


2 \ ho E — fo 


2 

V r = 


%Pl-Afoa 2 )PtJ A , 8fo 


( 


4«2 r 4 \0*E-p 0 
3J 2 fofo 


[PSE-fo 


1 3 (/3i — /3 3 2 /4q 2 ) v: 

*(p h-p° 


X 




(17) 


2 a 2 r 2 [PS E-fo 

In the case where /3 2 = Afo a 2 (the equivalent of 7 = a 2 for the quartic metric), the effective potential simplifies to: 

1 1 o 

3J 2 fo 


2p?v 2 r = 


PS E-fo) - 


0‘S r 2 [PS E - Po 


= -v eS , 


(18) 


and a relation between the potentials po and a 2 for v C i rc -C 1 may be determined by examining the circular orbits of as it was done 

for the quartic potential. If it is further assumed that a 2 is constant, the time potential Po satisfies a differential equation similar to that of the 
quartic potential /3 © showing that Po oc Mag log(r/ro) for an isolated mass. 


4 LIGHT DEFLECTION 


Our program of searching for static isotropic Finsler spacetimes for MOND based on the Tully-Fisher relation has so far led only to the 
determination of the metric’s time potential in the weak field approximation. Space potentials may be constrained from light deflection as 
the process is directly influenced by the space curvature of the metric (see below). For instance, in the Schwarzschild metric of an isolated 
mass in a Riemann spacetime, light deflection in the weak field limit is due in equal measure to time and space curvatures. In what follows, 
we examine light delfection of MOND’s Finsler metrics for small accelerations in the weak field approximation. The small acceleration 
approximation does not mean that photons are not relativistic. Similarly, the weak field approximation means only that the perturbation that 
deflects the free relativistic motion of photons is small (Schneider et al. 1992). 

For the quartic metric, light travels on the null surface given by F q = 0 = Pvf — (7 — a 2 )v 4 — ( av 2 + r 2 n|) 2 at constant energy 
E = Pv^ and constant angular momentum J = r 2 v < f,(av 2 + r 2 v 2 ). Its orbit 0(r) may be obtained from the three relations by eliminating 
v ^ as: 

J 4 /3(7</> ,_4 + r 4 + 2 ar 2 (j)'~ 2 ) 3 — E 4 r a (a<f>'~ 2 + r 2 ) 4 , (19) 

where the prime denotes d r . If the light surface is quadratic (7 = a 2 ) then the deflection equation is identical to that of Riemann spacetimes. 
Otherwise (7 ^ q2 )> < F~ 2 * s the solution of a polynomial of degree 6. For impact parameters b = J/E > roe -1 / 4Ma °, we may use the 
closest approach distance, r m defined as dr/d/6(r m ) = 0 and the notation /3 m = /3(r m ) to cast the deflection equation as: 


Prtni 7</>'" 4 + r 4 + 2ar 2 (j>'~ 2 ) 3 = P m r 8 (q ^'" 2 + r 2 ) 4 . 


( 20 ) 


We determine light deflection with the assumption that the potentials a and 7 remain constant as mass vanishes and spacetime may be 
described by a flat Lorentz metric in that limit. Writing a = 1 + a, P = 1 + /3, 7 = I + 7 , /3 m = 1 + pm where a,P, 7 <?C 1, and expanding 
( 120b to first order with respect to the perturbation, the deflection equation reads: 


(P-Pr, 


Tmu{u 2 — l ) 1 / 2 4 r m (u 2 — l ) 3 / 2 


(37 — 4a)u 2 + 3(2a — 7 ) 
4 r m u 3 (u 2 — l ) 1 / 2 ’ 


( 21 ) 


where it = r/r m . The three terms correspond respectively to rectilinear motion in flat empty space, deflection caused by time curvature and 
deflection caused by space curvature. 

Applied to the quartic Finsler spacetime of an isolated mass, the deflection produced by the time potential p = AagM log(r/ro) is 
A (f> = TvMao/2 similar to that from a singular isothermal sphere: linear with respect to mass and independent of the characteristic scale ro. 
Time curvature of the Finsler metric therefore cannot account for the larger observed deflections. F or co mp arison, interpretations of MOND 
based on gravity modification provide a deflection scaling Ac j> oc (Mao) 1 / 2 1 Milgrom 2013 : Mortlock & Turned 20011) that explains some 
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but not all of the observed mass discrepancy in galaxy clusters. For Finsler metrics, space curvature may account for a larger deflection in 
particular because of its possible dependence on the distance scale ro that in turn may depend on M and possibly a o. For instance, if we 
arbitrarily set a = c a Mao log(r/ro) and 7 = c^Mao log(r/ro), deflection is given by: 

A (j> = nMa ° [4 + - 2c a + (2c 7 - 4c a /3) log(6/2r 0 )]. (22) 

o 

In terms of the two distance scales of the system r u = a^ 1 , a constant related to the radius of the observable Universe, and r s = 2 M, 
the Schwarzschild radius, deflection A <j> oc r s log(&/2ro)r“ 1 for large impact parameters. The proportionality constant must be positive if 
deflection is to be convergent; this in turn, depends on the definition of r 0 whether it is larger or smaller than the intermediate acceleration 
radius r e = (r u r s / 2) 1//2 where the Newtonian gravitational acceleration is comparable to ao. Substituting galaxy-galaxy lensing typical 
parameters show that the arbitrary form of a and 7 that yields HD does not produce sufficient lensing whether ro is taken of the order of 
r s or r u . Obviously, this indicates only that the arbitrary choice of the space potentials is inadequate. Light deflection may constrain but not 
determine the space potentials as the effect essentially probes the integral potentials and not the corresponding acceleration through such 
integral equations as ED- We note that the deflection equation for the rational metric 03 is quite involved but its linearization is similar to 
ED and the time potential /3o yields a constant deflection for an isolated mass under the assumptions specified at the end of Section 3.2. The 
source of larger deflection for the rational metrics is also space curvature. 


5 DISCUSSION 


In our search for an interpretation of MOND as a modification of inertia at small accelerations, our starting point was Milgrom's original 
modification of Newton’s second law for circular motion. We have been able to show that implementing the Tully-Fisher relation in that 
case is simple but comes at a price. Although the gravitational potential remains linear with respect to mass, it does not satisfy the Poisson 
equation in the weak field limit as one would think from Milgrom’s original scaling by equating v 4 r~ 2 ag 1 and Mr -2 . The new gravitational 
equation 0 differs from a Poisson equation in that the divergence is not applied to the gradient of the potential j3 but to the product of d r /3 
and r u r _1 where r„ = a^ 1 is proportional to the radius of the observable Universe. There is also a difference in the meaning of the Finsler 
and Riemann gravitational time potentials in that the former is not equivalent to the latter but to the latter's square as may be seen on the 
quartic metric with a quadratic light surface 0. Although unexpected, the necessary modification of gravity preserves the linear character of 
the gravitational potential which may prove a precious advantage with respect to modified gravity theories. Moreover, covariant theories of 
modified gravity MOND require various auxiliary fields in addition to the physical metric in order to implement the Tully-Fisher relation. In 
this sense, our modification of inertia and its corresponding linear modification of gravity are minimal. 

Another element that appears in our analysis of asymptotic motion is the characteristic distance ro. Whereas it does not influence weak 
field massive motion, it is likely to contribute to light deflection through its possible dependence on mass and the acceleration scale ao- A 
covariant MOND formulation possesses two independent distance scales: r u , a multiple of the radius of the observable Universe, and r s , the 
Schwarzschild radius. It is not clear yet how the distance ro depends on these scales but its form should be related either to the transition 
from the small to the large acceleration domains or perhaps, more interestingly, to the presence of a coupling of local dynamics to large 
scale structure evolution through the acceleration scale ao - In this regard, our approximation of static and asymptotically flat metrics (section 
4) that implements the Tully-Fisher relation locally for weak field massive motion may not be adequate for light propagation. Although 
complicated by the absence of a full modified-inertia gravitational equation, it would be desirable to determine Finsler metrics for MOND 
equivalen t to the McVittie-type solutions of the Einstein gravitational equation for a mass embedded in an expanding Universe (M cYittie 
1933: lNandraetalj2012h , 

The appearance of the distance s cale rp has another implication regarding the possible scale invariance from which MOND may emanate 
in the domain of small accelerations. Milgrom ' 20091 20141 ) showed that the basic phenomenological MOND laws are natural consequences 
of the assumption that the equations of motion of purely gravitational systems are invariant under the change of ( t , r) into (At, Ar) when the 
acceleration scale ao tends to infinity. Our formulation that is based on implementing the Tully-Fisher relation indicates that the converse 
may not be true. In effect, the time potential we derived /3(r) = 1 + AMao log(r/ro) (for the quartic metric) enters the geodesic equations 
directly through the Finsler function so that the equations of motion are not scale invariant. The distance r 0 appears in various places: it 
is related to the distance at which the metric becomes singular (/3 = 0), although that distance is located far inside the Newtonian domain 
where the Finsler metric in its MOND form does not apply (lacking space and time cross terms). The distance ro also enters the expressions 
of the energy, angular momentum, and (proper) period. Light deflection may also depend on ro and in this respect our formulation departs 
significantly from modified-gravity MOND theories where asymptotic light deflection is constant akin to that of dark matter isothermal 
sphere models albeit with a different mass dependence. Constant asymptotic light deflection is sometimes considered a tenet of MOND 
phenomenology but its observational validity is questionable when applied to galaxy cluster encounters. We surmise that the appearance 
of the distance scale ro is ubiquitous in the context of MOND as whatever type of interpretation is considered (gravity modified or inertia 
modified) when formulated in a metric setting is bound to have an asymptotic metric that requires ro because of the logarithmic behavior of 
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the gravitational potentials. Once the scale r o enters a metric, the geodesic equations of motion obtained from the proper length action may 
no longer be scale invariant and light deflection will depend on the characteristic scale roQ 

The simplest Finsler metrics that are compatible with the Tully-Fisher relation make use of polynomials of degree 4 with respect to local 
displacements. The classical physical setting where such expressions appear is surface area measurements albeit it as the norms of differential 
2-forms. As the corresponding tensors are antisymmetric, we may look for symmetric 4-tensors in the context of area-metric spacetimes. 
An area metric is an invertible partly antisymmetric covariant tensor field, G« x^, satisfying: G K a m „ = G^^x, G^x^ — —Gand 
G riXah G ahlw = 2(5*5* — 5*5*). It gives rise to the totally symmetric Fresnel tensor field (J m „ k a defined by 


GfivnX — 


|detG| 

24~ 


^ijkl ^mnpq 


(~*ijm (/x G | fcn | k G \ ) Ipq 5 


(23) 


where e^^x is the totally antisymmetric tensor and G mi/ „a = G^^x — G[^ K x] is the cyclic part of the area-metric tensor. In a Riemann 
spacetime with a metric tensor q av , th e naturally induced area-metric tensor G M „ K a = g^nQvX — g^xg^n, and the Fresnel tensor simplifies to 
GuvnX = Q(,„,Q k \^ l Punzi et al.l2007 ). The Fresnel te nsor first appeared in the st udy of electromagnetism in anisotropic media l Tammll924 , 


925l : Rubilar 2002) and parti cularly of birefringe nce i Obukhov & Rubilail2002 ) as well as in the study of area-metric spacetimes that do not 


possess a prior metric tensor ' Punzi et afll2009bl) . Light propagation in the geometrical optics limit is shown to oc cur (in anisotropic medi a 
and area-metric spacetimes) on the null quartic surface defined by G tl ''**k l ikvknkx where is the wave covector I Hehl & Obukho\ll20Q3 >. 
The Fresnel tensor in an area-metric spacetime is therefore a possible candidate for constructing a Finsler distance function for MOND. The 
study of massive motion_in area-metric _spacetimes i s tric kier than light propagation as the fundamental objects on such spacetimes are not 
particles but strings i Schuller & Wohlfarthl 20061 ). Punzi et al. (2009a ) devised an ad hoc method for defining a localized massive particle 
by averaging the motion a string fluid and showed that such motion is described by the geodesics of a quartic metric defined by the Fresnel 
tensor. Understanding massive motion in area-metric spacetime is a first step towards determining the physical conditions that give rise to 
MOND-like inertia and how it is related to large scale structure evolution. 
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